A Representation for the Anyon Integral Function 



M. Aslam Chaudhry 1 , Amer Iqbal 2 , Asghar Qadir 3 

1,3 Department of Mathematical Sciences 
King Fahd University of Petroleum and Minerals 
Dhahran 31261, Saudi Arabia 

2 School for Mathematical Sciences 
Government College University 
Lahore, 54600, Pakistan 

2 Department of Mathematics 
University of Washington 
Seattle, WA, 98195, U.S.A. 

3 Centre for Advanced Mathematics and Physics 
National University of Sciences and Technology 
Campus of College of Electrical & Mechanical Engineering Rawalpindi, Pakistan 

Abstract 

The Fermi-Dirac and Bose-Einstein particles satisfy corresponding statistical dis- 
tributions. In the phenomena of charge fractionalization and the fractional quantum 
Hall effect it is found that particles behave as if they are neither fermions nor bosons. 
Such particles are called anyons. The integral functions for bosons and fermions are 
available in the literature. However, there is no anyon integral function available. In 
this note we propose a pair of functions that interpolate, in some sense, between the 
gamma function and the zeta function, which we call "gamma-zeta" functions. It is 
pointed out that this pair of functions very naturally provides a representation of the 
anyon integral function. 
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1 Introduction 



There are two aspects to this note: (1) the presentation of a new special function; and (2) 
its physical application. It is necessary to briefly review the background for both aspects, so 
that readers familiar with one area can follow the motivation and reasoning for the other. 
In this section, we will first recall the special function aspect and then go on to explain the 
physical aspect. 

Most of the so-called "special functions" arise as solutions of commonly occurring first 
or second order linear differential equations. However, two of the most significant special 
functions do not come from differential equations, namely the gamma function and the family 
of zeta functions. In fact both arose from the study of numbers. (For a discussion of the 
relevant special functions, see for example |TJ.) Considering their relevance and utility, and 
the fact that they are related in the above sense, it seemed worth while to try to put them 
together so that they appear as limiting cases of either one of a pair of functions, which 
we call the gamma-zeta function |2j- It turns out that in doing so we find that we can 
"interpolate" between the Fermi-Dirac and Bose-Einstein integral functions in a way that 
follows the physical needs for doing so, as will be explained. 

The gamma function was defined by Euler to extend the domain of the factorial function 
from the integers and then the complex numbers. It has the integral representation 



which is singular for negative integer values of a, but is well defined everywhere else. The 
zeta function, defined over the reals, was first studied by Euler who also gave a product 
formula for it known as the Euler product, 



Riemann while studying the distribution of prime numbers rediscovered the zeta function 
and studied it as a function defined over complex numbers. He also established the amazing 
link between the complex zeroes of zeta function and the distribution of prime numbers. The 
Riemann zeta function, as it is now known, has the integral representation 
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The function has the analytic continuation |3] 
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is the loop integral. The contour C consists of the real axis from oo to p (0 < p < 2ir), 
the circle \z\ = p, and the axis from p to oo. The integral 1(a) is uniformly convergent 
in any finite region of the complex plane and so defines an integral function of a. The 
representations (J3J) and (JHJ) provide an analytic continuation of ((a) over the whole complex 
plane. Since 1(1) = 2i\i and J(n) = (n = 2,3,...), the poles of the gamma function 
T(l — a) at a = 2, 3, 4, . . . cancel the zeros of 1(a). It follows from (JSJ) that ((a) has only a 
simple pole, due to T(l — a), at a = 1 with residue 1. 



2 Bose-Einstein and Fermi-Dirac Particles 

In quantum mechanics particles (for example see jlj) are represented by wave functions that 
satisfy the Schrodinger equation. Making the theory relativistic requires that one shift over 
to quantum field theory (for example see Following the Dirac quantization procedure 
leads to the Klein-Gordon equation, which has problems of interpretation due to the fact 
that it has a second derivative relative to the time appearing in it. To avoid these problems 
Dirac "took the square root of the Klein-Gordon equation" to obtain a first order equation in 
space and time variables. It turned out that the Klein-Gordon equation represents spin-less 
particles while the Dirac equation applies to particles with spin. The theory further requires 
that the spin go up in half-integer multiples of the quantity, h, defined to be h/2iT, where h 
is Planck's constant. The wave function for identical particles with integer spin is symmetric 
while for half-integer spin it is anti-symmetric under the exchange of particles. On account 
of this property, particles with half-integer spin can not have the same quantum numbers, 
while those with integer spin can. The former are said to be subject to the Pauli exclusion 
principle. 

Statistical mechanics (for example see [6.) had been developed to deal with ensembles of 
classical particles using the Maxwell distribution: 

f(e) = e^l\ (1) 

where r stands for thermal energy and is hence related to the temperature, e stands for 
the kinetic energy of the given particle and p for the chemical potential. The / gives the 
probability of the particle having the given kinetic energy. It turned out that this distribution 
did not apply precisely to the quantum particles but could be used as an approximation for 
a mixture of the two types. Also, it was a good approximation for the ensemble at high 
temperatures. For lower temperatures of systems of particles of one type other distributions 
were required. For the integer spin particles it was propounded by Bose and is called the 
Bose-Einstein distribution: 

Mc) = e(c J /T _ r (2) 

while for the particles with half-integer spin it was propounded by Fermi and is called the 
Fermi-Dirac distribution and is given by: 

If(s) = -, 1 -. (3) 

2 



Particles of the former type are called bosons and of the latter type fermions. 

The cumulative probabilities for bosons and fermions are given by the Bose-Einstein and 
the Fermi-Dirac integral functions: 



i r°° t q 



■= tttt—^ I -zzr^rdt (? > 0); (4) 



1 C°° t q 

FJx) := r / dt (q > -1). (5) 



3 Anyons 

It appeared that all elementary particles would belong to one of these classes. However, more 
recently, in the phenomenon of the fractional quantum Hall effect [7j it was found that under 
certain conditions electron, a fermion, can behave as if it is made of yet more fundamental 
particles with spin a fractional multiple of \h. Such particles were dubbed "anyons" jH] 
(For a review see j^]). Particles such as anyons with spin/statistics interpolating between 
bosons and fermions can only exits in two dimensions. There are both topological and 
group theoretic reasons for their non-existence in higher dimensions. As mentioned in the 
last section, the wave function of bosons is symmetric under the exchange of two particles 
whereas the wave function of the fermions is antisymmetric. This has to do with the fact 
that in three and higher dimensions the symmetry group is the permutation group and it 
has only two one dimensional representations. The trivial representation corresponds to the 
bosons and the non-trivial representation gives fermions. In two dimensions, however, the 
situation is much more interesting. The symmetry group is larger than the permutation 
group, it is the braid group 4 . The braid group has a one dimensional representation for 
every real number v. Thus a two particle wave function under exchange of two particles 
behaves in the following way: 

{(+1) xjj(a,b), Bosons, 
(—l)ip(a,b), Fermions, (1) 
( — l) u ip(a, b), Anyons. 

There is no functional representation for anyons available in the literature corresponding 
to that for bosons or fermions. The gamma-zeta functions connect across the two in that 
for real u they give an interpolation between them. As such, they should be considered as 
candidates for representing the anyon integral function. Note that $ and \& are sufficiently 
directly related that either (or a linear combination of the two) could be used. The choice 
would be guided by examining how convenient the function is in the resulting formulae. 

It was found ^U] that there are phenomena in which it appeared that there were fractional 
charges appearing. This was regarded as exciting because it seemed that these may be the 
"quarks" that had been proposed in high energy physics ^T]. However, the theory did not 

4 Braid group, -B/v, is generated by transpositions {4>i, • • • ; 4>n} such that <j>i<j>j = <t>j4>i f° r \i ~ j\ > 1 and 
4>i<pi+i4>i = 4>i+i ( f > i4>i+i- If we further impose the relation <f>f — \ then we get the permutation group. It is 
the absence of this last relation in the braid group which allows the possibility of anyons 
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allow free quarks jT2] and it was found that the fractions did not necessarily correspond to 
those required for quarks. This phenomenon was later explained as arising from the statistics 
not being well defined, namely that the particles behave like neither fermions nor bosons but 
something in between. The explanation proposed can be understood in terms of an example. 
Consider fermions of spin 1/2 (in units of the Planck action h) in an infinite linear array. 
The axial component of the spin can then either be +1/2 or -1/2. The energy is minimized 
by having them alternate in the axial component of the spin. Had they been bosons there 
would be invariance under shift by a lattice length. For fermions we would have to reverse 
the sign on the shift. This can be achieved by shifting and then "rotating the particle" 
through an angle it, i.e. changing the phase of the particle wave function. This way we can 
have the fermion behave like a boson in some sense. This is relevant for our purposes. More 
generally the transformation may be through some other fraction of 2n. 

We should be able to follow the same procedure for the integral functions, of shifting the 
phase by 7r, or more generally, by some other fraction of 2n. Our proposed new "gamma- 
zeta" function has the Bose-Einstein and Fermi-Dirac integral functions as limits. As such, 
it should be able to achieve just this. As will be seen, the procedure adopted corresponds to 
the physical argument just mentioned. 



4 The Gamma- Zeta Functions 

Our pair of new functions, which we call gamma-zeta functions, is obtained by using the 
Mellin and Weyl transforms. The Mellin transform of a function <p(t), if it exists, is defined 
by US! 

POO 

$ M (a) := M[<p; a] = / * a- V(t)<ft (a = a + i6). (1) 

Jo 

It is inverted by 

i pc+ioo 

cp(t) = — / * u {z)t-'dz. (2) 

/7T£ J c-ioo 



A convolution-like product of two functions ip and ip, <p o tp, can be defined by \L£ 

POO 

(cp o ip)(t) := / cp{xt)ip(x)dx. (3) 



Note that the operation "o" is not commutative and 

M[<poil>;a]=$ M (a)V M (l-a), (4) 

where 

$ M (a) :=M[<p;a] (5) 

and 

tf M (a) :=M[V;4 (6) 
It is to be noted that if (p G -£q oc [0, oo) is such that 
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then $m(«) exists for o\ < a < Oi- In particular, if ip is continuous on [0, oo) and has rapid 
decay at infinity, the Mellin transform (4.1) will converge absolutely for a > 0. 

Let Xm(cti, cr 2 ) be the space of all ip G L\ oc \fi,oo) such that the corresponding integral 
(4.1) converges uniformly and absolutely in the strip <j\ < a < oi. Then, for each ip G 
Xm{&i, &2), ®m(oi) is analytic in the interior of the strip <j\ < a < cr 2 , [IS] - 

Similarly the Weyl transform of ip G -X"jy[0, 00) is given by [Tlj 

$ w ( a; x) := W- a [ip\{x) = ^^M[<p(x + t);a} (a = a + zr, a > 0), (8) 

where 

:=<p(x), (9) 

and 

$ w (-a;:r) := D m ($ l4 /(a; x)) (a = a + it, o > 0), (10) 
where m is the smallest integer greater than a = Re(a) and 

D n :=(-l) n £; (n = 0,l,2,...). (H) 

In particular we have 

$ w (-n; *) = (-l)>W(s) (n = 0, 1, 2, . . .). (12) 

For our purposes we need to consider some sufficiently well-behaved functions. The class 
£(<7 ) of "good" functions is defined as follows: A function ip G S(cr ) if 

(P.l) <peX M (0,a ), 

(P.2) $v^(a,x) is analytic in x in the region Re(x) > for < o < a . 

It is to be noted that the class £(o"o) is nonempty. To see this consider any polynomial P(t). 
If a > 0, the function P(t)e~ at is a member of the class £(<7o)Vcxo > 0. Also, as <3>iy(0, a;) is 
analytic in the region Re(x) > 0, each of ip G £(<7o) is a C°°[0, 00) function as well. 
Let us now define 

e l — 1 



and 

ibJt) : = 

e* + l 



^(*):=^-TT. (^>0,t>0), (14) 



and consider their Weyl transforms: 

<S> u (a;x):=W- a [if u }(x) (15) 

and 

V v {a;x):=W- a [^]{x). (16) 
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We call the functions $ u (a;x) and *Sf v (a;x) the gamma-zeta functions. Note that 



3& 



$v{at]x) = t^t r / t a - l Vv {t + x)dt (17) 

-ra /"oo ^a— lg— ft 

-dt 



r(a) 7 e*+* - 1 

-ft 



r(«) y e f -e 1 

Similarly, we have 



^• x) = ^Wl VT^ At - (18) 

From (J 17)) and (|18|) we find that the gamma-zeta function $^(0;; x) is related to the general 
Hurwitz-Lerch zeta function ^3] &{ z , s, v) via 

$„(a; x) = e' {v+1)x ^(e~ x , a,u + l), (19) 

which leads to the series representation 



n + u + l) a 



n=0 

(v ^ -1, -2, -3, . . . ; Re(x) > 0; x = 0, a > 1). (20) 
Similarly, we have the relation 

x) = e- {u+1)x <5>{-e- x , a,u+l) (21) 
that leads to the series representation 



4r„(a; x) = e"^ V / iJ 6 , 



^ / -1, -2, -3, . . . ; Re(x) > 0; x = 0, a > 0). (22) 



Theorem 1 TTie functions $> u (a;x) and ^ u (a; x) are related via 

*„(a; x) = (-l) (i/+1) $,(a; x + nr). (23) 

Proof. Replace x by x + m in (}2T)j) and compare with (}2"2*|) to obtain the result. ■ 
Another useful relation, also proved elsewhere [2] is: 

Theorem 2 The functions $> u (a;x) and $? u (a;x) are related via 

tf 2 „(a; x) = $ 2i >; x) - 2 1 - Q $,(a; 2x). (24) 

6 



Proof: 



n>0 

2 e -(2k+l)x 



{2k + 2 + 2u) a 
2 i-a e -xJ2 ' 



k>0 

-2kx 



n+ 1 + v) a 

k>0 v ' 

= 2 1 - a e- x $ v {a-2x)e {u+1)2x 

Note that our gamma- zeta functions are dual to each other in the sense that the above 
relation can be easily inverted, so that each is similarly related to the other. This fact is of 
special relevance for us. 

5 The Anyon Integral Functions 

From the integral representations (j!7j) of § u (a;x) and (2.4) of the Bose-Einstein integral 
function B a (x), we find that 

B a _x(-x) = $ (a;x), (1) 

leading to the fact that the gamma-zeta function x) is a natural extension of the 

Bose-Einstein integral function. Similarly, from (fTH|) and (2.5), we find that 

=# (a;x), (2) 

which shows that the second gamma-zeta function ty u (a;x) naturally extends the Fermi- 
Dirac integral function. 

Putting v = and replacing x by — x and a by a + 1 in (|24j) we find that the Fermi-Dirac 
and Bose-Einstein integrals are related via 

F Q (x) = B a (x)-2 1 - a B a (2x). (3) 

The relation (5.3) does not seem to have been realized earlier. 

Using the fractional Weyl transform W~ a on <& u (f3;t) we obtain 



1 



oc 



$„(a + (3- x) = — - I t a ~^M t + x). (4) 

r(a) Jo 

Putting v = above we find the interesting integral representation 

1 f 00 

B a+ p_ l {x) = -— 1?*- x Bp-. x (x-t)dt (5) 
1 l a J Jo 

for the Bose-Einstein integral. Similarly, we obtain 

1 f°° 

J- i«J Jo 



(6) 
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Putting a = 1 in (5.5) and (5.6) we obtain 

/•oo 

B p {x) = / B p ^(x - t)dt, (7) 
Jo 



and 



poo 

F p (x) = / F . x [x-t)dt. 
Jo 



These representations for the Bose- Einstein and Fermi-Dirac integral functions may prove 
useful. 

We can interpolate between the Bose-Einstein and Fermi-Dirac integral functions using 
gamma-zeta function. To see this consider G u (a; x), 

G u (a, x) = a(v)® v (a; x) + 6(i/)* I/ _i(a; x) . (9) 

If we appropriately chose a(v), b(u) satisfying the conditions a(0) = 1 — b(0) = 1 and a(l) = 
6(1) — 1 = then G u (az;x) interpolates between the Bose-Einstein and Fermi-Dirac integral 
functions. 



6 Concluding Remarks 

The family of zeta functions including the Riemann, Hurwitz, Epstein, Lerch, Selberg and 
their generalizations, constantly find new applications in different areas of mathematics 
(number theory, analysis, numerical methods, etc.) and physics (quantum field theory, 
string theory, cosmology, etc.) A useful generalization of the family is expected to have wide 
applications in all these areas as well. On the mathematical side, the gamma-zeta functions 
$„(a; x) and ^ v {a\ x) discussed in this paper provide a unified approach to the study of the 
zeta family that is remarkably simple. This is discussed in more detail in a separate paper 2j. 
For our present purposes of the physical applications, the first pair of gamma-zeta functions 
also provide a representation of an anyon integral function. There was no such representation 
available in the literature. They also led to a new relation between the Bose-Einstein and 
Fermi-Dirac integral functions that reflect the nature of the two functions. The "mixing" 
of fermions and bosons displayed in this relationship demonstrates how the fractional spin 
behaviour arises. 

One might have hoped that the gamma-zeta functions could be applied to supersymmetry 
16 , as they "unify" bosons and fermions in some sense. Despite the "duality" manifested 
by them, they cannot be so used. The reason is that supersymmetry does not interpolate 
between the two types of particles. They remain quantized with spin whole or half integer 
with nothing in between. The utility of the gamma-zeta functions is that they do provide 
such an interpolation. This lack of intermediate particles is in sharp contradistinction to 
the situation for grand unified theories ^7j, in which quarks and leptons are unified. In 
this case there are intermediate particles, called lepto- quarks that are neither quarks nor 
leptons but both. Such particles are required for spontaneous symmetry breaking of the 
unified symmetry ^Hj- This highlights the problem with supersymmetry, that it does not 
have a spontaneous symmetry breaking mechanism and has to use "soft symmetry breaking" . 
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Again, one sees a strong relationship between the physical and mathematical requirements 
in the context of our gamma-zeta functions. 

It is expected that the operational properties of these functions can and will be exploited 
further in solving difficult problems in mathematics and physics. 
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